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Abstract 

In the first part, expressions are given for the sign of the topolog- 



o 

f—^ . ical angle that is acquired upon making a loop around a degeneracy 

("conical intersection") point of two molecular energy surfaces. The 
expressions involve the partial derivatives (with respect to the nuclear 
coordinates) of the matrix elements of the coupling Hamiltonian. Ex- 
amples are given of a few studied cases, such as of excited states that 
C3 have topological angles with a sign opposite to those in the ground 

states. 

In the second part, the two dimensional (or two parameter) situ- 
ation that characterizes a conical intersection (ci) between potential 
surfaces in a polyatomic molecule is constructed as a limiting case 
of the three dimensional Dirac-monopole situation. For an electron 
occupying a twofold state, we obtain both the " magnetic- field" (or 
curl-field) and the tensorial (or Yang-Mills-) field (which is the sum of 
a curl and of a vector- product term). These pseudo- fields represent 
the reaction of the electron on the nuclear motion via the nonadia- 
batic coupling terms (NACTs). We find that both fields are aligned 
with the orthogonal, (so called) seam directions of the ci and are zero 
everywhere outside the seam, but they differ as regards the flux that 
they produce. In a two-state situation, the fields are representation 
dependent and the values of, e.g., the fluxes depend on the state that 
the electron occupies. The angular dependence of the NACTs and the 
fields calculated from a general linearly coupled model agrees with re- 
cently computed results for C-xB. [A.M. Mebel, M. Baer and S.H. Lin, 



J.Chcm. Phys. 115 3673(2001)]. An effective-Hamiltonian formalism 
is proposed for experimentally observing and distinguishing between 
the different fields. 



1 Introduction 

In 1975 Longuet-Higgins showed that if there is a sign change in the wave 
function upon performing a closed loop in the parameter space, then a point 
of degeneracy is enclosed by the loop [1]. Trivially, a sign change amounts to 
a phase angle change of ±(2iV + l)ir, where TV is an integer or zero and is, in 
fact, zero when the circulation is about a degeneracy point which is a single 
conical intersection (ci). In the following section of this work we consider 
the sign of the added phase angle change, first in a general manner, then for 
a special ("the complex") representation, including several examples that 
have been discussed in the literature ([2] - [10]). The sign associated with 
any single ci will be of added interest in cases that there are several ci's in 
the system. We shall focus attention on such multiple ci situations. The 
next section has been motivated by some recent publications that connect 
the electron-nucleus interaction with the three-dimensional Dirac-monopole 
formalism [11] and with a pseudo-magnetic field [12] due to the nonadiabatic 
coupling. In the latter paper an angular dependence of the field was postu- 
lated. We now derive and specify the unique form of this field for the most 
general model having a linear form of electron-nuclear interaction near the 
conical intersection. A further pseudo-field (the Yang-Mills tensorial field) 
is also introduced in the present molecular context and is evaluated. 

2 Signs of Geometrical Phases around Conical In- 
tersections 

2.1 Cartesian, real representation 

Let X and Y denote the coordinates of the parameter-plane in which the 
circling is being performed and let (Xq,Yq) denote the location of a ci. We 
shall consider the case where the potential surfaces of only two states of the 
system cross at this ci. In the vicinity of the ci we are entitled to disregard 
possible interactions with all other states and to write the potential energy 
V(X, Y) as a 2 x 2 matrix that has the following form. 
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Figure 1: Coordinate systems for the point (X,Y) or (q,(f>) that circles 
around the point of conical intersection [the dot located at (Xq,Yq) or 
(qo,(f>o) ]. The circling is with a radius A and circling angle a. 
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In the absence of a magnetic field the matrix components can be taken 
to be real. A scalar term has been omitted, since this would not affect our 
considerations. At a ci we have 



A(X ,Yo) = 0, B(X ,Yo) = 



(2) 



but the partial derivatives of either A or B with respect to the coordinates 
(to be designated by subscripts X and Y) cannot both be zero at (Xo,Yq). 
(Otherwise, the intersection is touching, not conical.) The mixing angle in 
the eigenstates of the potential matrix equation (1) is given by 



\ arctan l^X) 

-arctan-—- 

2 ysC 



X )B X (X ,Y ) + (Y- Y )B y (X ,Y ) 
X )A X (X , Y ) + {Y- Y )A Y (X , Y ) 



(3) 



the second line being valid close to the intersection point. Use has been 
made of equation (2) . We now introduce the radius of circling A and the 
circling angle a shown in the figure, for the case that Ax is non-zero. 

(If Ax = 0, then Ay must be non zero and we use an angle oriented 
by 2^ from that shown in the drawing.) The introduced quantities are 
implicitly defined by 



X — Xq = A cos a, Y — Yq = A sin a 



(4) 



We now note the important fact that for a ci, in its vicinity, the acquired 
phase is monotonic with the circling angle and it is therefore permissible to 
obtain the sign of the phase by looking merely at the limits 

a -► (5) 

or, equivalently, 

cos a— ► 1, sina^a (6) 

We thus find for the phase angle, after some elementary simplification, 

6 « - arctan(— + a -p ) ( 7 ) 

with all derivatives to be evaluated at the ci. Since the inverse tangent is 
an increasing function of it argument, as the circling angle a increases, the 
mixing angle 6 will increase or decrease depending on whether 

A X B Y - B X A Y > or < (8) 

(We recall that Ax was assumed to be non-zero). 

Example 1: 

A pair of conical intersections, e.g., (a) in some bent molecules, as H2S 
of Civ symmetry, two ci's between the lowest excited states 1 Ai and 1 Bi, 
(b) for C 2 H in C s symmetry the 3 2 A' and A 2 A' states [13, 14] 

We now write out the Hamiltonian representing the coupling between 
the nuclei, whose coordinates are denoted by (X, Y) and the two electronic 
states whose surfaces are intersecting conically. For the case of a pair of ci 
symmetrically situated in the X, K-plane the situation can be described by 
the following matrix elements: 

A(X,Y) = (X 2 -l), B(X,Y) = Y (9) 

We have so shifted and scaled the nuclear coordinates, that the ci are 
located at X = ±1,Y = 0. The representation symbols of Herzberg [15] 
have been used, for which X transforms in Civ as &2, Y as B\ and Z as 
A\. We verify that at both ci the partial derivative Ax is non vanishing. 
Furthermore, since 

A X B Y - B X A Y = 2X (10) 

has opposite signs for the two ci, we immediately find that the topological 
angles for the two ci have opposite signs. Therefore, circling around the two 



ci gives a total accumulated phase of zero, rather than 2tt, as might have 
been expected. This was indeed the result found computationally for large 
radius circling in some low lying 2 A' states in C 2 H [16]. Thus we conclude 
that the above model is suitable for this molecule in the vicinity of the ci. 

Example 2: 

In C 2v symmetry the degeneracy between the 2 B 2 and 2 Ai states in 
Al — H 2 [3] . Here symmetry arguments forbid the presence of an off diagonal 
term of the form shown above in equation (9) . Noting further that the ci 
are in a plane orthogonal to the plane formed by the three atoms, we use 
(X, Z) instead of (X, Y) and write the off diagonal term as B(X, Z). Terms 
that are made allowed have the form 

A(X, Z) = (X 2 - 1), B(X, Z) = XZ (11) 

The reason that this off-diagonal term is allowed, is that X transforms as 
B 2 and Z as A±. The two ci are located at X = ±1, Z = 0. We now find 
for the Jacobian 

A X B Z - B X A Z = 2X 2 (12) 

and this gives the same sign for either ci. The acquired phase is therefore 
the same and circling around both yields 2tt. 

The results of both Examples 1 and 2 have been previously obtained by 
us, using a different, graphical-algebraic procedure, the "continuous phase 
tracing method" [17]. 

2.2 Complex representation 

This representation is commonly used, e.g. for the doublet electronic states 
in the E<S>e Jahn- Teller problem [18, 2]. The cartesian representation (X, Y) 
of the nuclear displacement coordinates is conveniently replaced by the cylin- 
drical polar coordinates (q,(ft), according to X = qcos(ft,Y = q sin (ft. It has 
been shown [19] that when the vibronic coupling (between nuclear and elec- 
tronic degrees of freedom) is taken to any arbitrary order, then the matrix 
in equation (1) (expressing this coupling) can be put in a complex form, 
such that its off-diagonal part has the form 

V l2 {qA) = Kqe~^[l+q~ 2 £ q* m Q m+ (q)e Zmi <t>+ £ q^Q m ^q)e^ mi( t>\ 

m=l... m=l... 

(13) 
Here K is a constant and Q m + and Q m - are polynomials in q 2 with real 
coefficients that depend on the physical system (or of the model used to 



represent it) and whose leading term will be q°. (An alternative expression in 
a real representation is found in [20].) Normally, for stable physical systems, 
it is expected that with increasing m, Q m + and Q m - will both numerically 
decrease and so will, in each polynomial, the coefficients of successively 
higher powers of q 2 . 

There will also be diagonal, scalar and pseudo-scalar terms. The scalar 
term simply shifts the energies of the two states by an equal amount and 
can be ignored. The pseudo-scalar term has the form 

J2 q^ m P m {q)sin?,m^ (14) 

m=l... 

where the polynomials P m are defined in a similar manner to the Q's, above. 
This term arises when the system is not time-reversal invariant, such as in 
the presence of a magnetic field. We shall not consider this situation. By 
consequence, the level separation is simply 2\Vi2(q, <f>)\ and the mixing angle 
is 

9 = — arctan (15) 

2 ReV 12 

Now, conically intersecting degeneracies of the two (complex) electronic 
states can occur (those ci are not the only ones possible) at two types of 
points: At q = 0, coming from the initial, q- factor in equation (13) and 
at trigonally located degeneracies (qo,<f>o) from the square-bracket factor in 
Vi2(q, (/>), such that one has either 

0o = 0, 2tt/3, 4tt/3 (16) 



or 



vr/3, vr, 5tt/3 (17) 



To find the sign of the phase acquired upon circling around any trigonal 
ci, we use, as above, the circling radius A and the circling angle a shown in 
the figure. The polar coordinates, expressed in terms of these are, correct 
to the first order in A, 

q w qo + A cos(a — </>o), ~ <fio H sin(a - 4> ) (18) 

We shall now find it convenient to rewrite equation (13) in the form 

V 12 {q,(t>) = e-^U 12 {qA) (19) 



where Ui2(q,<j>) is defined by 

U 12 {q,<j>)=Kqe-^-^[l + q~ 2 £ q 3m Q m+ (q)e 3m ^ 

m=l... 

+ E <? 3m Qm- (g)e- 3 ™ ] (20) 

m=l... 

When this quantity is expanded about the ci, one obtains by equation (18) 

At t a £\t j 

U\2{q,4>) ~ U 12 {qo,(/>o)+Acos(a-(i)o)(^—) qo ,4, +—sin(a-(f)o)(-^—) q0j( t >0 

oq go o<p 

(21) 
The first term on the right vanishes at a point of degeneracy. Real and 
imaginary parts of Ui2(q,<j>) come from the middle and the last terms, re- 
spectively. Therefore, by equation (15) , 

j / dUi2 \ 

9 = -i<p + o arctanf - — ,t,°' ° tan(a - (f>o) ) (22) 

The sign of the term linear in the angle a — <po will determine the sign of 
the topological phase. 

Example 3: Let us choose, as in an earlier work [10], a quartic approxi- 
mation for the off-diagonal matrix element 

V 12 (q, <j>) = Kqe-^1 - fiqe M * + A^e" 3 ^] (23) 

We first consider the case in which \x and A have the same sign. Then, in 
addition to the zero at the origin, this expression can have up to three sets 
of trigonally positioned zeros, each of which gives rise to three ci. This is the 
maximum number of zeros, but since a negative or complex q^ is physically 
not admissible, there may be less than 3 x 3 ci. 

To determine the sign of the phases, one searches for the coefficient in 
equation (22) of the linear term in a — </>o as follows: At q = 0, the coefficient 
is a constant: — g- This gives, upon full circling, the usual topological phase 
of — it. At a trigonal point, one finds for the prefactor of a — 4>o 

3^90 + 3Agg 
Hqo - 3Agg 

With the parameters A and /j, having the same sign, the above ratio is positive 
or negative, depending on whether 

\ — > or < qo (the radial position of the ci ) (25) 

V 3A 



When equality holds, there is a double root and the intersection is no longer 
conical. 

With the choices of [x = 0.3, A = 0.003 (used for illustration in [10]), the 
following roots are found: 

One root at q = 

Further, the following trigonal roots: 

go = 3.95,0o = 0, 2vr/3, 4vr/3 
go = 7.42,0o = O,27r/3,4vr/3 
q = 11.37, (p = vr, vr/3, 5vr/3 (26) 

Substituting the parameter values in the left hand side of the inequalities in 
equation (25) , one finds for the square root the value of 5.77. This makes the 
signs of the phase of consecutive sets of ci come out positive, negative and 
negative (in the order of increasing radius), resulting in accumulated phases 
of -7r (around the origin only), then 2ir, then -it and ultimately, at large 
radius circling, -4n. All of these values were confirmed independently by the 
"continuous phase tracing method" [10] and by numerical computation of 
the phase accumulated along the angular coordinate. In the model of [2], and 
of others ([3], [4] -[6]) identical to that one, the off diagonal matrix element 
stops with the quadratic term. In the language of our model, exhibited in 
equation (23) , this means that the coefficient A of the cubic term is zero. 
In virtue of the criterion expressed in equation (25) this implies a positive 
phase around each of the three trigonal ci, or for large-radius circling, when 
the central ci is also circled, a net phase change of +2ir. This agrees with 

the results of the cited works. 

i 

The case of different signs in a and A and also I — j-l < 1 yields three 
new type of roots, as follows: 

% = \f^ (27) 

V A 

required to be real and positive, with angles 4>o at the intersection that are 

1 




3 arccosh/— /(2 M )] (28) 

in which the argument of arccos is required to be magnitude-wise less than 
one, as well as two further angles oriented with respect to the previous by 



2n/3 and 4tt/3. The roots can be checked upon substitution into equation 
(23) • 

The signs of the phase angles can be obtained from equation (15) and 
equation (21) . After some manipulations one obtains an expression, which 
reads, correct to the first power of tan(a — <fio), 

11 3 

6 = — arctan[- tan(30o) + -(tan(3</>o)) 2 £an(a — 4>o)] (29) 

Since arctan is an increasing function of its argument, and noting the neg- 
ative sign in the expression, we find that circling around these ci gives a 
negative phase angle. This last type of ci, which is also trigonally located 
but angularly shifted from the positions given in equation (16) or equation 
(17) , has not been noted before, neither experimentally, nor theoretically 
or calculationally. 

2.3 Phases in excited states 

(a) A three dimensional parameter space. 

The following Hamiltonian was used by Berry [7] for elucidation of the 
geometric phase in a three-dimensional parameter space represented by 
R = (X, Y, Z). The same Hamiltonian formed the basis of discussion of 
a magnetic monopole in real space [8]. 



(30) 

This can also be written in terms of Pauli spin matrices, in the form 

H=^(Xa x + Ya y + Za z ) (31) 

The (adiabatic) eigenstates of the Hamiltonian can be obtained in terms of 
the spherical coordinate representation of R = (R, 6, 4>) as 

^={ cos***/* ) (32) 

for the lower state, and 

^ = I sinleW ) (33) 




for the upper (excited) state. The separation of the states is A = R. It 
was shown, in e.g. [9], that the geometric angle 7c arising from describing a 
closed contour C in the R- (or parameter) - space can be expressed in terms 
of the expectation value of gradient of the Hamiltonian 



A n H = -(±a x + Ya y + Za z ) (34) 



namely 

r MA R ff|vo 
icm = J — ^2 dS ( 35 ) 

where the integral is taken over a signed area S enclosed by the contour C. 
Evaluating the expectation values and carrying out the surface integration 
for a contour at constant 9 (namely, over a region enclosing a spherical cap), 
one obtains the following geometric phases for the two eigenstates. 

7cC0*) = -(1 -cos#)7T (36) 

for the lower (ground) state, and its opposite 

7c(V«) = (l-cos0)7r (37) 

for the upper (excited) state. One sees that, for a contour around the large 
circle, for which 9 = ^n, the phase-factors, namely e l7c , are the same (—1) 
for the two states, though not the phase- angles, whereas for other values of 
the angle 9 not even the phase factors are the same. In principle, phase factor 
differences can be established experimentally by interference measurements. 

(b) Adiabatic (slow) time development in a doublet. 

We consider the time-dependent Schrodinger equation written as 

i^A = H {x ,tMx,t) (38) 

(in which t is time, x denotes all particle coordinates, H(x, t) is a real time 
dependent Hamiltonian and % = 1). As is well known, the presence of i 
in the equation causes the solution ^f(x,t) to be complex-valued. We now 
consider a special case for the Hamiltonians in which two externally imposed 
independent sinusoidal perturbation interact linearly with the electron. Fur- 
ther, the electronic Hilbert space is confined to a degenerate electronic dou- 
blet, represented by two vectors: ( J ) and [ ° ). In this representation the 
Hamiltonian can be written as a 2x2 matrix, for which we have chosen the 



10 



following form, whose physical origin has been described in various works 

(e.g., [21] - [23]) 

H(t) = G/2(-™W S[n ^\) (39) 

Here u is the angular frequency of the two external disturbances, taken to 
be the same for both. The eigenvalues of equation (39) are — ^ and ^. We 
take G > 0, so that the former is that ground state energy. A method of 
solution was outlined in [23] and the amplitude of the vector ( J ] in the 
ground state was given there by the expression 

x f (t) = cos(Kt) cos(ujt/2) + (u/2K) sin(Kt) sm(cot/2) 

+i(G/2K)sin.(Kt)cos(ut/2) (40) 

with 



K = 0.5VG 2 + lo 2 ^G/2 (41) 

the latter approximation being justified in the adiabatic, slow motion limit 
defined by 

G/uj > 1 (42) 

In the present work we do not repeat the method of solution in [23] , but use 
it to present the results also for the amplitude of the vector ( ° ) for the 

ground state to be denoted by xf(*)> as weu as the corresponding amplitudes 
X\{t) and xlW i n the excited state. The two states are differentiated by the 
initial conditions, for t = 0: 

Xf(0) = l,xi(0) = (43) 

for the ground state and 

Xf(0) = 0,xl(0) = l (44) 

for the excited state. We recall the well known result that in the adiabatic 
limit the state maintains its ground or excited state character throughout 
the motion [24]. After some algebra one obtains the following result, which 
encapsulates the four amplitudes in one single formula for arbitrary initial 
values Xi(0) and X2(0): 



Xf + A (t) - ( _\ )^{e^-^[xi(0)(^ + f + |)(l( t )/ie 



iu)t\ 
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+ z X2 ( )(if-| + |)(l(;)/ 2 e-*)] 

e - i{K -iu,)t [ximK _ | + |)(i( + )/ 2 e--*) 

- i X2 ( )(tf + | + |)(l( ; )/ ie -- *)]} (45) 

where the upper and lower symbols in the parentheses have to be taken 
consistently with the choice of the components. The formula, as presented, 
is valid generally, for arbitrary boundary conditions. For the special choice 
of either the ground or the excited state the initial conditions, as expressed 
in equation (43) or in equation (44) , have to be substituted in the above 
formula. The auxiliary functions f\ and ji are given by 

K -i- Q — t ±L K — Q- — *iL 

nT 2 T 2 2^2 

In the adiabatic limit, when G/u 3> 1 and therefore K = ^G 1 + to 2 ~ 
G/2 correct to the first order in %, the expression simplifies to take the 
form: 



+ <X2(0)(^)(i( ; )/ 2 0] 

( t ) e-( G -^[xi(0)(^)(l( + )f 2 e-**) 



%(0)(1 + — )(1(; j/ie— *)]} (47) 

The crucial element in the establishment of the proper sign of the topo- 
logical phase is the recognition that the prefactors f± and / 2 of the circular 
functions e ±la; ' in equation (46) are magnitude-wise less than unity. There- 
fore, upon circling a full period, none of the round brackets containing these 
in the above expression can cause a sign change. Thus the sign of the phase 
is dictated by whether the first or the second term in equation (47) stays 
finite (non-vanishing) in the adiabatic limit % — ► 0. It is now evident from 
the initial conditions in equation (43) that for the ground state it is the first 
term that stays finite, and therefore the topological phase is — n, whereas for 
the excited state, for which the initial conditions are equation (44) , it is the 
second term that survives and yields the topological phase of +ir. Note also, 
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that the first exponential factors e arise from the dynamical phase and 
are irrelevant to the topological phase in the adiabatic limit. For their effect 
upon the topological phase in the not fully adiabatic case, we refer to [25]. 
We further find from the expression that, in either the ground or the excited 
state, both components 1 and 2 (diferentiated in the above expressions by 
round and square brackets, respectively) have the same topological phase. 

Thus, in summary, we find that the signs of the topological phases are 
identical for both components in a given state, but are opposite in the ground 
and excited states. It is possible to interpret this result from a wider per- 
spective, as follows: 

The time-dependent ground state of a bound system is characterized by 
a lower boundedness of the energy. This leads to the analyticity of the wave 
function (more precisely of its logarithm) in the lower half of the complex 
t-plane [26, 27]. The phase is then obtained from the modulus by a Kramers- 
Kronig type relation, which involves integration along a contour consisting 
of the real t-axis and a large semicircle in the lower half of the t-plane 
[22, 23, 25]. In contrast, the upper partner in a doublet has an energy upper 
boundedness. This leads to an analyticity in the upper half t-plane, to a 
closing of contours (in the Kramers-Kronig relation) in the upper half plane 
and to a topological phase of the opposite sign. 

In a general situation that involves in the Born-Oppenheimer superpo- 
sition more than two electronic states, one still expects a negatively signed 
topological phase in the adiabatic ground state, but in general one cannot 
predict the sign in an intermediate-energy adiabatic state. 

3 Pseudo-fields Arising from Degeneracies 

Denoting by \i > a set of electronic adiabatic eigenfunctions and by V the 
gradient operator conjugate to the nuclear coordinates, we introduce the 
non-adiabatic coupling terms (NACTs) as 

Tij =< *|V|j > (48) 

It was noted some time ago that the NACTs can be incorporated in the 
nuclear part of the Schrddinger equationas a vector potential A [28]- [32]. 
The question of a possible magnetic field H has been considered in [12] . 
This field is not due to any source external to the molecule, but rather arises 
from the mutual coupling between the electrons and the nuclei. To mark this 
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point we call H a pseudo-field. Formally, it is associated with A through 

H = curl A (49) 

In the present work we derive the magnetic field H, for the case where 
the electron- nucleus coupling is expected to have an especial importance, 
namely, for two (diabatic) electronic states becoming degenerate at a sin- 
gle point in a two dimensional (2D) parameter space. In the context of a 
polyatomic molecule (which is our subject of reference), this 2D space is the 
plane created by a pair of nuclear displacement modes, to be designated X 
and Y. The degeneracy of the electronic states gives rise to a ci, convention- 
ally regarded as the origin X = Y = 0. The NACTs have a pole at the origin 
and, excluding this point of singularity but still staying in the neighborhood 
of the ci , the curl of the NACTs vanishes [34, 35]. Farther away from the 
ci , the curl of the NACTs is non- vanishing, provided some ci due to other 
adiabatic states are found in the region. These two facts (the "vanishing of 
the curl" around the ci and the existence of a pole at the origin) cause the 
issue of the magnetic field to be problematic. The approach in this work 
is to obtain the field associated with the ci by a limiting procedure, which 
yields the field unambiguously. 

The magnetic field that we derive is tied to the fact that in the neighbor- 
hood of the ci we have an electronic sub-manifold of at least two dimensions. 
(The actual electronic Hilbert space is, of course, of higher dimensionality, 
but this higher dimensionality can be ignored near the ci , whereas a one 
dimesional space cannot give rise to a ci, by definition.) However, as will 
be verified in the sequel, the NACT-induced magnetic field will vary with 
the kind of superposition one makes of the electronic states within the sub- 
manifold. In the following section we first obtain the magnetic field in the 
representation of the adiabatic states. These are the appropriate choices 
under slow variation of nuclear coordinates, since then the adiabatic states 
are eigen-states which are uncoupled from the rest of the manifold. Subse- 
quently, we shall also consider a magnetic field in another (non-adiabatic) 
representation. The representation dependence of the magnetic field is con- 
nected to the non-Abelian nature of the situation. (" Non-Abelian" means 
the involvement of more than one electronic state and the presence of non 
commuting matrices in the Hamiltonian. The need to specify the represen- 
tation does not arise in Abelian (commuting) systems, which is the usual 
background for classical electromagnetism. There one defines a unique mag- 
netic field without regard to the electronic state the system is in. The mag- 
netic field is, trivially, "unique" for a single electronic state. However, it is 
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such also for a many state situation, such that the Hamiltonian contains no 
non-commuting matrices.) 

Another objective of this work is to obtain, in addition to the magnetic 
field, the tensorial Yang- Mills field. By its original conception [36] and its 
numerous later developments and applications (e.g., [7] and [38]), the Yang- 
Mills field constitutes a residual interaction within the manifold (in our 
context, between the electronic states). From a historical perspective, there 
is a long and distinguished line of papers in which a magnetic field has been 
obtained for a related problem: a pair of states that become degenerate at a 
point in a three dimensional parameter space. Notable among these works 
are Dirac's derivations of the quantized monopole- field [39], Wu and Yang's 
matched vector potential [40] and Berry's formulae for the magnetic field 
(e.g., [7], [41], as are also other papers reprinted in Shapere and Wilczek's 
volume [42].) 

3.1 The adiabatic states 

We start with the following Hamiltonian (taken over from Berry's ground 
laying paper [7] with one essential modification) 

This is written in a diabatic representation of the electronic states (that 
are such that the state does not involve, at least locally, nuclear coordinates). 
A Hamiltonian matrix written out in a different representation, employed, 
e.g., in [12] but formally different from that in equation (50) is considered 
later in this paper. The two components involved in the formalism are 
denoted by |+ >,|— >, and the nuclear mode coordinates by (X,Y, Z)(= 
R). The parameter b in equation (50) is to be noted. The magnetic 
monopole treatment (e.g., in [7, 40]) has 6 = 1 and "spherical symmetry". 
We shall obtain the NACT and a magnetic field for a general b ^ 1, and 
then proceed to the 6^0 limit, which is the situation of interest to us, 
namely the two parameter problem. (The parameter Z varies along the so- 
called seam-direction of the degeneracies. In the limit of b = 0, Z does not 
enter the dynamics. Physically, bZ represents a magnetic field aligned with 
the quantization direction of the two states in equation (50) . We write it 
as ±bZ and then let 6^0, rather than let Z — ► 0, first to maintain the 
analogy with [7], and also since a coordinate Z cannot be made to vanish in 
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a Hamiltonian.) The matrix in equation (50) has eigenvalues ±R b and the 
corresponding adiabatic states take the form (for b > 0): 

|1 >= e -^/ 2 cos(^/2)|+ > +e^ 2 sin(fl 6 /2)|- > (51) 

|2 >= - e -^/ 2 S in(0 6 /2)|+ > +e l(p/2 cos(0 6 /2)|- > (52) 

where 



R b = ^JX* + Y 2 + {bZf = y/q* + (bZy (53) 

tan (j) = Y/X (54) 

tan 9 b = q/bZ (55) 

having introduced the radial coordinate q in the cylindrical polar description 
(q, cp) through 

q 2 = X 2 + Y 2 (56) 

The NACTs have already been defined as elements of the matrix (de- 
noted by T h ) of the gradient operator V conjugate to the variables (X, Y, Z). 
Thus t is both a vector and a matrix. Its elements in the adiabatic-state 
representation, equation (51) - equation (52) , are: 

t\ x = < 1|V|1 >=--icos6 b V<p (57) 

= -r h 22 (58) 

Tl 6 2 = ^sin0 6 V<£-^V0 6 (59) 

T% 1 = -ism6 b V<f>+-Ve b (60) 

[This is a departure from the usual definition in the chemistry litera- 
ture of NACTs, which are real off-diagonal and anti-symmetric matrices 
(or matrix elements). However, we shall find that our diagonal NACT are 
anti-hermitean, just as are the usual NACTs. (Trivially, i times the NACTs 
are hermitean. The usual definitions of the NACTs can be regained by re- 
defining, say, the second adiabatic state states to be i times the ones used 
here. Our choice makes the "off-diagonal" magnetic field, which can be seen 
below in (16) and (18), come out real.) Particle physicists use the term 
connectivity (A), which also is a matrix- vector, related to r by 

A = ir (61) 
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(as in [37]). They further employ the terms curvature (or the Yang- Mills 
field intensity tensor), defined by 

F = VAA-iAAA = H-iAAA (62) 

and the Berry-phase or line integral over a closed contour C 

A • dR (63) 

c 

(The conditions, in a molecular context, for the vanishing of F have been 
described in [34].) 

The magnetic field H arising from the NACTs is 

H = iV A r (64) 

We first calculate this in the representation of the adiabatic states shown in 
equation (51) -equation (52) . 

H$i = -nl 2 = iVAT b u (65) 

= - V A (cos 6 b V<t>) (66) 

= sin0 6 (V0 6 A'vV) + -COs6> b (VA'vV) (67) 

and 

H 6 12 = H b 21 = iVAr b 2 (68) 

= — V A (sin 9 b V(/)) (69) 

= — cos6 b (V6 b AV</») - -sin6> 6 (VAV(/>) (70) 

having put V A V9 b = by elementary vector algebra. The (apparently) 
similar quantity V A V0 needs a different treatment, because of the singular 
nature of <fi on the seam line q = 0. This will be done below. 
We next calculate the Yang-Mills (tensorial) field, given by 

F = i(VAr + TAr)=H + irAr (71) 

also in the representation of the adiabatic states 

F b n = HJ 1+ 4at^ (72) 
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1 -> -.1 - 

= -V A{cosd b V(j)) + -sme b {V6 b AV(l)) (73) 

= --sm0 b (V0 6 A V</>) + -cos0 b (VA V0) 

+ isin0 6 (W 6 AV0) (74) 

= -cos6» fe (VAV0) (75) 

= "F^ 2 (76) 



and similarly 



¥\ 2 = HU^Ar^^Ar^) (77) 

= -~VA(sin0 b V^) + -cos0 6 (V0 6 AV</>) (78) 

= -cos0 5 (V0 6 A V0) - -sin0 6 (V A V</>) 

- -cos0 b (V0 b AV(/>) (79) 

= --sin0 6 (VAV<?!>) (80) 

= F^ (81) 

In the above expression we call attention to the cancellation between the 
terms in the curl and in the vector product (the first and the third terms in 
equation (74) and equation (79) ) 

The derivative quantities V<?5> in the NACTs (equation (58) and equation 
(59) ) and V A V</> in the field intensities (equation (75) and equation (80) 
) do not depend on b. They have been considered in [32], section III.C in 
a somewhat different context, namely as arising from a "pure- gauge" phase 
factor that multiplies the whole adiabatic state. Their conclusion is that they 
give rise to a pseudo-magnetic field that is zero everywhere except along the 
"curve of intersection [the seam] of the two potential surfaces ...where it has 
a delta- function singularity". We shall use this result in the formulae that 
immediately follow, but note that quantities that depend on b are new to 
the present work and so are some results that do not vanish in the limit of 
6^0. (A formal justification of the result of [32] involves the extension of 
Stokes' theorem to singular integrands and will be given elsewhere [43].) The 
(^-derivatives will now be given in two coordinate systems: 1) a cylindrical 
coordinate system (q, <p, Z) with unit vectors depicted by bold and hatted 
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symbols and 2) in the Cartesian space (X,Y,Z), for which we quote the 
results of [32], using the notation of (i, j, k) for the unit orthogonal vectors. 

V0 = 4>/q [= (Xj - Yi)/(X 2 + Y 2 )} (82) 

V A V0 = Z5{q)/q [= 2nS(X)5(Y)k] (83) 

S(q) is the Dirac delta function. This result is in accord with [34]. 

It thus immediately follows that the tensorial fields shown in equation 
(76) and equation (81) are zero, except on the seam, irrespective of the 
value of bj 

3.2 The Non Adiabatic Coupling Terms (NACTs) 

The sin 6b and cos 6b prefactors of these derivatives depend on the parameter 
b, as the notation indicates, and the results are true for any, general value 
of b. However, we shall examine what happens in the 6^0 limit, namely 
as the 2-parameter (or 2-dimensional) problem is reached. We recall that 
the angle 6b has been defined in equation (55) . From this definition, in the 
limit b — ► 0, we have, for q ^ 0, 

limsin^(= q/R b ) - 1 (sin(|) = cos(|) = -L) (84) 



Clearly, also 



lim cos 6 b {= bZ/Rb) -► ir5{q)R b -+ n5{q)q (85) 

6^0 



To obtain the last result we have used the expression for the delta func- 
tion 

nS(q) = lim — 5 ^- (86) 

and the definition of Rb in equation (53) . We add two further results to be 
used in the sequel. The first is 

V6 b = -(bq/R 2 b )Z + (bZ/R 2 )q [= -(bq/R 2 b )k + (bZ/qR 2 )(Xi + Yj)] (87) 

and from this, using equation (82) , 

V0 6 AV0 = (b/R 2 b )ci+(bZ/qR 2 b )Z (88) 

[= (b/qR 2 b)(Xi + Yj) + (bZ/qR 2 b)k] (89) 

= (bZ/RlM/Z + Z/q) (90) 
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Then, from equation (57) and equation (82) , 

A = - % -{bZ/qR h )^ (91) 

[= - l -( bZ /R b )(Xj-Yi)/(X 2 +Y 2 )] (92) 

= -A (93) 

h i 1 - bZq , q Z . , , 

i. i 1 - bZq , q Z N , N 

74 = 2** + ^7"^ (%) 

From equation (85) we see immediately that (because of the delta func- 
tion factor arising from bZ/Rb) the diagonal NACTs are zero outside the 
seam line q = as b — ► . However, on the seam line they are non-zero. 
This is a new result, which has been obtained by going to the limit 6^0 
after calculation of the derivatives. 

If we now calculate the Berry-phase or the line integral shown in equation 
(63) taken with a circular contour encircling the seam line at any finite 
distance (q > 0) from it, we obtain 

An ■dK = TTcos9 b (96) 

c 

In the 2-parameter limit, as b — > 0, this will tend to zero (from equation 

(85) ). By familiar arguments (based on Stokes' theorem that equate the 

contour integral with a surface integral of the curl), the flux across any 

finite part of the parameters plane of the "diagonal" magnetic field in the 

adiabatic representation vanishes, notwithstanding the fact that the vector 

potential is not zero inside the (infinitely thin) q = "solenoid". (We shall 

presently check this result by actually calculating the surface integral of the 

"diagonal" magnetic field.) In the off-diagonal NACT given in equation (59) 

or equation (60) , the second term can be seen [from its form in equation 

(87) ] to give S(q)q as b —> 0, whose line integral around a finite circle 

q > clearly vanishes. However, the first term in the off-diagonal NACTs 

survives the 6 — *• limit even for q > 0. Thus we have a genuine, finite 

vector potential entering through the off-diagonal part. The resulting line 

integral 

/ Ai2-dR = -7rsin6» 6 (97) 

Jc 

is responsible for the Berry phase of -it. 
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3.3 The magnetic field in the adiabatic representation 

As already noted above, we can check the previous line-inegral results by 
the equivalent method of evaluating the flux of the magnetic field H across 
an XY- plane. H can be calculated from equation (67) and from equation 
(70) . For the diagonal part we obtain 

H 6 n = [(bZ/2qR b )5(q)-(bZ/2R 2 )}Z-(bq/2R 2 )q, (98) 
[= (n(bZ/R b )6(X)S(Y) - (bZ/2R 3 b ))k 

-(b/2R 2 )(Xi + Yj)} (99) 

= -H 6 22 (100) 

The second expression (in square brackets) is written in the Cartesian frame. 
In the 6^0 limit this field vanishes outside the seam line q = 0. In the 
general case, for which b > the last term is oriented along the cylindrical 
radius vector (q) and only the Z component of the field contributes to the 
flux across the (q, <fi) [or the (X, Y)]-plane. Corresponding to the result for 
a circular contour at q shown in equation (96) which vanishes in the 6^0 
limit (as seen above), one obtains after some manipulation of the surface 
integral 



lim(M agnetic flux) n = n dq5(q) — irbZ / dq—pr (101) 
6^0 ' Jo Jo R h 

= (102) 

= — (Magnetic flux)22 (103) 

Thus the two methods (the line integral and the surface integration) give 
the same result. 

The off-diagonal magnetic field is from equation (70) 

[= (-7r(q/R b )5(X)5(Y) - (bZ) 2 /2qRl)k (105) 

-{b 2 Z/2qRl){Xi + Yi)\ (106) 

= H*i (107) 

In the limit 6^0, this field also vanishes outside q=0. To evaluate the 
flux, care must be taken to go to the limit b — ► only after the integrations 
are performed. At the end, one again obtains equation (97) . 
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3.4 The Yang-Mills fields and collected results 

These are defined in equation (71) , are shown for the adiabatic representa- 
tion in equation (76) - equation (81) ) and have the following forms: 



and 



Fix 


= (bZ/2qR b )5(q)Z 


(108) 




[= (n(bZ/R b )5(X)6(Y))k] 


(109) 




— ~*22 


(110) 


F fe 12 


2 R b 

[= -7rq/R b )5(X)5(Y)k] 


(111) 




= F b 21 


(112) 



These are much simpler than the magnetic (purely curl) field expressions. 
Their fluxes are given below, in Table 1, where we collect all results in the 
2-parameter 6^0 limit. 



Table 1. Summary for conical intersections in the adiabatic representation 
in the two-dimensional, 6^0 limit (designated by sub- or super- script.) 
Derivatives: V0 = 4>/q, V8 = 7r%)q, V6» A V0 = nS(q)(q/Z + Z/q) 
NACT: rft = -r 2 ° 2 = -(i/2)7rS(q)4>;^ 2 = r& = ^/q - i<5(g)q 
Flux or line integrals: 

(Magnetic/ lux) n = — (Magnetic f lux)22 = 
(Magnetic/ lux) 12 = (Magnetic/ lux) 2 i = — vr 
(YM/lux)ll = -(YM/lux)22 = vr 
(YM/lux)l2 = (YMflux)2l = 



As already noted, the most remarkable results in Table 1 are the non- 
zero values of the NACTs upon the seam q = and the circumstance that 
the magnetic and Yang-Mills fluxes appear in an opposite manner. These 
are the legacy of our using the 6^0 limiting procedure, rather than starting 
with a Hamiltonian in which 6 = 0, the Z-coordinate is absent and the seam 
q = is a line of singularity. The observational significance of the fields or 
fluxes will be discussed later. 
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3.5 Circulating representation 

In the foregoing, we have calculated the field and the flux on the assumption 
that the electron is in an adiabatic state. For a slowly changing perturba- 
tion, such a state is a stationary state of the system. The diagonal elements 
of the tensorial quantities that we have calculated refer to this situation. 
In principle, the electron may be excited to a state different from an adia- 
batic one, e.g., to a linear combination of two adiabatic states with constant 
coefficients. The superposition of adiabatic states will persist, since under 
quasi-stationary conditions each adiabatic state will develop in an indepen- 
dent fashion. The fields pertaining to such a situation can be obtained 
through use of the non-diagonal elements of the field (or flux), calculated 
in this work. However, the magnetic field can also react on the electronic 
motion and this reaction might cause changes in the adiabatic state. This 
effect is formulated in the last subsection. 

We now calculate the fields in a different representation, which we name 
"circulating representation" and denote by circularly shaped bras and kets, 
|+),|— ). This is the representation introduced by Baer [12] so as to ob- 
tain NACTs that are diagonal (pseudoscalar) and purely imaginary. (It is 
suggested that in a non-reactive scattering, the |+), |— ) states resemble the 
circulating-state situation.) 

The following transformation, acting on the adiabatic states in equation 
(51) , equation (52) , generates the circulating representation: 



|+) = [|1>+|2>]/v / 2 (113) 

|-) = [|1>-|2>]/V2 (114) 



In the limit 6^0 (where sin -^ = cos -y = l/\/2), one has 

|+) = e 4|-> (115) 

|-) = e - i l|+> (116) 

where, we recall, |— > and |+ > are diabatic coordinate independent elec- 
tronic states. The complex exponential prefactor is the reason for the name 
"circulating representation". [When b is non-zero, the circulating states are 
coordinate dependent superpositions of the diabatic states.] We list the 
results in this representation, obtained after using some algebra: 
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Table 2. Results in the circulating representation with the arrow denot- 
ing the 6^0 limit: 
NACTs: _ 

b _. J4_ i± 

'++ 2R b 2q 

T b — i<L 

T — ~ 2R b 

b _ ibZ4> 1 bZ a i 1 bqZ 
T -+ ~ 2qR b 2rJ^^ 21% 

b _ ibZ4> | 1 bZ a 1 bqZ 
T +- - 2qR b "T 2Rj*l 2"r2" 

Fluxes or line integrals, for b — > 0: 

(Magnetic/ lux) ++ = — (Magnetic/ lux) — ► — tt 

(Magnetic/ lux) _| = (Magnetic/ lux) (_ — v 

(YM/lux) ++ = -(YMflux)— -» 
(YM,/7ux) + _ = (yM/Iux)_ + -» 7T 



We note that for b nonzero the NACTs have also off-diagonal elements 
and that on the seam q = these elements exist even in the two-dimensional 
limit 6^0. Moreover, these elements have components other than tangen- 
tial (along 4>). As such, they do not contribute to the line integral in equation 
(63) . At finite distances from the seam q > 0, our results agree with those 
in [12]. 

3.6 An alternative formalism 

Our results could also have been obtained, had we proceeded differently, 
namely by making the X and Z the "active" variables and Y the fictitious 
one. To achieve this one puts the b factor (whose limit is ultimately taken) 
in front of Y in the Berry Hamitonian, our equation (50) , rather than before 
Z. This amounts to a representation of the diabatic electronic states that is 
obtained from the |+ >, |— > set used above, by a complex transformation. 
In the 6^0 limit the adiabatic states are the same in either formalism and 
so are the NACTs and the magnetic field. An historical interest is attached 
to the latter procedure in that Stone not only proposed this Hamiltonian 
[46], with the interpretation of b as a spin orbit coupling strength, but even 
considered the limit of b — > . However, he did not derive a magnetic field. 
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3.7 General linear coupling and derived quantities 

The most general form for the interaction Hamiltonian that is linear in the 
parameters (X, Y, Z) is the following: 

^ aX + if3Y -bZ j v ; 

This differs from equation (50) by the inequivalence between not only Z 
(for 6^1) and X, Y, but also between X and Y. One regains the for- 
malism of [47], upon using the substitutions b = 0, a — ► b, (3 — ► 1. In this 
two-dimensional (or two-parameter) case the intersection between the two 
adiabatic potential surfaces (the solution of equation (117) ) involves two 
inverted elliptic cones, rather than two inverted circular ones which one ob- 
tains for b = 0, a = (3. The method of the previous section can be easily 
carried over to this case and we shall only quote the modifications needed 
for the generalization. Our main purpose is to demonstrate explicitly the 
dependence of the magnetic and Yang-Mills fields on the angle <p defined in 
equation (54) . A dependence of this type has been predicted in [12]. We 
first define 

7 = a/j3 (118) 

so that 7 = 1,(3 = 1 takes us back to the previous sections. Otherwise, 
(when 7 ^ l)we have the following modifications, designated by placing an 
apostrophe over all symbols affected: 

The adiabatic states have now the following form: 

|1' >= e-^'/2 cos (^/ 2 )|+ > +e ^'/ 2 sm(6' b /2)\- > (119) 

|2' >= -e-^'/2 S in(0// 2 )|+ > + e i *7 2 cos(0£/2)|- > (120) 

where 



R' b = ^{aXy + {/3YY + {bZ)* = sje + {bZ? (121) 

tancj)' = j3Y/{aX) (122) 

tan9' b = q'/bZ (123) 

having introduced the modified radial cylindrical coordinate q' through 

q' 2 = {aX) 2 + {(3Y) 2 (124) 
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The NACTs are now 






'b - 
T ll — 


< i'|v|l' >=• 


- - cos 6 


= 


'b 

r 22 




r 12 


= -isin^V^' 


-i?« 



(125) 
(126) 

(127) 



r'l= l -isiud' h V<t>' + l -V6' h (128) 

The derivatives of the angles are 

^' = 7—7 2 7 n 2-7^/9 = ^# (129) 

1 + (7 — 1) cos z <p q z 

^ A ^' = 7-?-^T7-^T— Z = ^<K<z)Z (130) 

1 + (7 — lj cos z <p q q z 

W£ = -(bq'/R b 2 )Z + (&Z/i#)^ - -^(« 2 - /3 2 ) Sin2# (131) 

q 2q'R b 2 

and 

V^ A V<j>' = f^(q/Z + Z/g) (132) 

In equation (131) the last term is new. In the adiabatic representation 
equation (119) and equation (120) one finds the following NACTs 



" = ->'«*> 1 + (7 » -1) erf » * <133) 



T- 

z 

>b 



= -t& (134) 

and 

-12 = (2^ + («-/ 3 ) Bm2 ^)^-2^(g-^ (l35) 

-2i = (W b ~ (a ~ P)sm2 ^ ) ^ + ^R? { - q ~^ ) (136) 

In the two parameter limit 6^0 for non-zero q the only term remaining is 
r[l = t£ (137) 

w>^ (138) 

(139) 



L 6 

i 



2 1 + (7 2 - 1) cos 2 4> q 
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The off-diagonal coupling term obtained in [13] computationally for the 
molecule CiH near a conical intersection was subsequently fitted to the 
above expression in equation (139) [47]. The coupling term is character- 
ized (for 7 < l)by two (frequently quite dominant) peaks at (j) = 0(= 2ir) 
and 7T. However, on the seam q = there are additional non-zero terms 
in t ± 2, in r^, etc., as we have already noted. These are new results. One 
expects similar angular behavior from them, too. The line integral over the 
coupling coefficient, or the Berry phase, in equation (63) gives in the 6^0 
limit the values of for the diagonal terms and ±7r for the off diagonal 
terms, irrespective of the value of the ratio 7 = a/ (3. These values were 
confirmed numerically to a good approximation for the neighborhood of the 
intersection in [13] and discussed in [47]. 

Returning to general values of the parameters b, a, (3, we show now the 
fields (magnetic and Yang- Mills): 

H'A = a/3[-^(^|%)-^)Z-(6 (7 /2 J R; 3 )q] (140) 

2K b q K b 

= -H£ (141) 



H'A = -^S+&4-^4 (142) 



<?%) , {bZ? )± b 2 Z q 
2q'R' b 2q'R' b 3> 2R' b 3 q' 

U' 2 \ (143) 



In the limit 6^0, this field also vanishes outside q = 0. However, 
on the seam q = we obtain one of the interesting results of this section, 
namely the angular dependence of the magnetic field. Only the second term 
in equation (142) for H 12 survives in this limit and one obtains after a slight 
simplification 

HS = -- , f {q \ 9 Z (144) 

12 2 1 + ( 7 2 -l)cos 2 V ; 

An angular dependence of the seam magnetic field was postulated in [12]. 
The above relation gives its form within the general (elliptic) linearly de- 
pendent model. The total flux in the 6^0 limit is again the limiting form 
of equation (97) , namely — it. 

The Yang Mills tensorial fields are calculated as: 

, b Q bZq5{g) ^ 

11 = ^ 2q' 2 R' *> ' 

= -F£ (146) 
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whose integrated flux in the 6^0 limit is n for all values of a and (5 . 
Finally, 
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-<t* <-> 



?'b 



= F 2 6 ! (148) 

whose flux is for all finite values of b. Since in our limiting procedure the 
limit b — ► is taken at the end, we have a zero flux in the two-dimensional 
elliptic geometry, which is identical to the entry in Table 1 for the circular 
case. 

3.8 Interpretation of the fields 

Associated with a ci of two potential surfaces for a polyatomic molecule, 
there exists an analogue of a magnetic field, which affects the nuclear mo- 
tion through its presence in the nuclear Schrodinger equation. However, in 
addition to the magnetic field, there is an analogous, symmetry-based field, 
the Yang-Mills or tensorial field. By employing a limiting procedure we have 
shown that for a (multiple valued) adiabatic state both types of the field 
have delta-function-like, thin solenoidal forms. 

In the molecular context, the direction of the solenoid is defined by the 
ci, as follows: The ci defines a plane in the nuclear coordinate space; the di- 
rection of the field is along any arbitrary direction in the nuclear coordinate 
space (the "seam"), which is perpendicular to this plane. Though the two 
types of fields have similarities, they differ in their numerical values and the 
fluxes due to them "complement" each other (Tables 1 and 2). We have also 
found that the fluxes are quantized, meaning that the strength of the field 
and the flux associate with it do not depend on physical parameters, just 
as they don't for the magnetic monopole field of Dirac [39]. This result was 
obtained recently [12], where the "curl- field" was calculated by applying a 
complex-valued linear superposition of the adiabatic states. This superposi- 
tion is equivalent to that introduced here for a pair of circulating states (one 
clockwise and another anti-clockwise) around the seam. It also expresses the 
geometry of the Aharonov-Bohm effect, in which two currents circle around 
a screened solenoid in opposite senses. However, in that case the flux is not 
quantized, but depends on the magnetic field inside the solenoid. 

As noted above, the "curl-field" is a delta function along the seam, so 
that a particle circulating at a finite radius would be oblivious of this field. 
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A different situation could arise when the electron is excited into a general 
superposition state, i.e. one that is not a superposition of adiabatic states 
(or is a superposition thereof, but with coefficients depending on the co- 
ordinates.) The "magnetic field" would be completely different from those 
for the adiabatic states. Formally, the new, general superposition would 
be described by applying what is called a "non-local" gauge transformation 
[37, 38]. The effect of this is well known and is expressed by saying that the 
vector potentials A (or the NACTs) transform inhomogeneously and the 
tensorial field F does so homogeneously in a covariant way [37]. However, 
for a coordinate dependent superposition, the fields may be difficult observe 
(if at all possible), due to the off-diagonal matrix elements in the potential, 
which make such a state non-stationary 

3.9 Observational aspects through effective Hamiltonians 

Possible experimental consequences of gauge fields have been noted for elec- 
tron spin experiments with time- varying magnetic fields [48], in atoms with 
rotating electric fields [48], in collisions between atoms [30] and in further 
applications [42]. 

We now give a general formalism for the observational effects of the 
fields, one that holds the promise of differentiating between the magnetic 
and the Yang- Mills field. It is based on an effective or truncated Hamiltonian 
formalism, similar in many respects to the well known Spin-Hamiltonian de- 
scription ([49] -[51]). This concentrates on a small set of states (in the 
present context, the two-fold set in equation (51) and equation (52) ), and 
considers the effects of perturbations on these. The perturbations admix 
states from outside the small set. The Spin-Hamiltonian formalism shows 
a way to represent the effect of the full set within the small ("truncated") 
set, in such a manner that the excited electronic states are included only 
"virtually". Symmetry considerations determine the form of the truncated 
Hamiltonian. Clearly, since several independent expressions can be com- 
patible with the symmetry requirement, as e.g., by including higher order 
effects, the effective Hamiltonian will normally consist of several terms. The 
coefficients with which these terms enter will in general not be amenable to 
calculations, but are empirically determined. 

After these preliminary remarks we recall that the previous subsections 
treated electronic and nuclear degrees of freedom. The fields that we have 
derived were pseudo-vector quantities in the nuclear space and were func- 
tions of the nuclear variables (X, Y, Z) or (q, 9, <fi). Moreover, the fields were 



2<s 



elements of a matrix (or tensor) in the Hilbert space of the electronic set, 
which is not the case in ordinary electromagnetism. We shall write a generic 
field component (not differentiating for the moment between magnetic and 
Yang-Mills fields) as 

F^ n = F^ n (X,Y,Z) (149) 

where a is the vector-index in the nuclear coordinate space and (m, n) are 
labels for the electronic set. We shall also consider operators of two types. 
First, the operators 

op1,opf,etc. (150) 

in the nuclear vector-space that are functions of electronic variables. It is 
now elementary to construct terms for an effective Hamiltonian H e ffi such 
that satisfy the symmetry requirements: 

H ef f = \m> (CiF^ r < r\opi\n > +C 2 F% ip Fp r < r\opf\n > +...) < n\ 

(151) 
with possible additional terms to follow. Summation for repeated indexes is 
implied. The first operator opf is a pseudo-scalar and the second operator 
0P2 has also the appropriate transformation properties. (Thus it might be a 
direct product of two opf s.) The coefficients C±, C 2 , etc. are expected to be 
empirical parameters. 

As a second application we consider the extension of the electronic (or- 
bital) degrees of freedom, e.g., by inclusion of electronic spins. We label 
the set for the extended degree of freedom by \M > and operators that act 
on both the electronic orbital and spin degrees of freedom by Op\,Op\ , 
etc. The a, b indexes refer to the axes in nuclear space, as before. Then the 
effective Hamiltonian takes the following form: 



H eff = \m>\M > (CiiC < r\ < M\Opf\n > \N > 

mp-"- pr 



+C 2 F*F b pr <r\< M\Opf\n > \N > +...) < n\ < N\ (152) 



Evaluation of these Hamiltonians requires computing expectation values in 
a given nuclear state. Since H e ff is anticipated to be small compared to 
energies of the nuclear freedom, this is a legitimate procedure. The effec- 
tive Hamiltonian thus provides a way to include residual perturbations that 
couple to states outside the degenerate electronic doublet. We have seen 
above (e.g. in Tables 1 and 2) that the magnetic field and the Yang- Mills 
field differ markedly. It is therefore suggested that by experimentally test- 
ing the effective Hamiltonian (e.g., through its dependence on the nuclear 
vibrational levels) one could establish which field is effective. 
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4 Conclusion 

This paper explains the signs of topological phases obtained by various au- 
thors in several previous works. The sign depends on the derivatives of the 
coupling matrix and is shown by the inequalities in equation (8) for a carte- 
sian, real representation and in equation (22) for a complex representation. 
Although positively and negatively signed topological phases cannot be dis- 
tinguished, since they have the same phase factor, this is true only upon 
performing a complete loop. For loops shorter or longer that this, the sign 
of the phase change is observationally accessible. 

We then obtain for various models the state- or representation- depen- 
dent magnetic and Yang- Mills fields, which result from the Born - Oppen- 
heimer scheme for the coupling between nuclear and electronic degrees of 
freedom. An effective or truncated Hamiltonian, here suggested, provides a 
possible tool for experimental verification of the fields. 
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5 Figure Caption 

Figure 1. 

Coordinate systems for the point (X, Y) or (q, <j>) that circles around the 
point of conical intersection [the dot located at (Xq,Yq) or (qo,4>o) ]• The 
circling is with a radius A, and circling angle a. 
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